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Abstract 

We investigate the quantum versus classical dynamics of a microwave cavity-coupled-Cooper 
pair transistor (CPT) system, where an applied dc bias causes the system to self-oscillate via the 
ac Josephson effect. Varying the dc bias allows the self-oscillation frequency to be tuned. An 
unusual feature of the system design is that the dc bias does not significantly affect the high 
quality factor of the cavity mode to which the CPT predominantly couples. The CPT-cavity mode 
system has a mechanical analogue involving a driven coupled pendulum-oscillator system. The 
corresponding, nonlinear classical dynamical equations exhibit chaotic, as well as aperiodic motions 
depending on the initial conditions and the nature and strengths of the damping/noise forces. The 
quantum master equation exhibits such phenomena as dynamical tunnelling and the generation of 
nonclassical states from initial classical states. Obviating the need for an external ac-drive line, 
which typically is harder to noise filter than a dc bias line, the self-oscillating system described 
here has considerable promise for demonstrating macroscopic quantum dynamical behavior. 



* To appear in Fluctuating Nonlinear Oscillators, Edited by Mark Dykman (Oxford University Press). 
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I. INTRODUCTION 



The work presented in this chapter has its origins in a seemingly mundane microwave 
engineering question: is it possible to apply a dc voltage (or current) bias to the center con- 
ductor of a superconducting coplanar microwave cavity, without significantly affecting the 
quality factor of, say, the first and second microwave modes of the cavity? Our original mo- 
tivation behind this question was to devise a circuit quantum electrodynamics (QED) based 

n n n 

scheme [1[ that can generate and detect quantum states of a mechanical resonator [2|, |3[, 
where the dc bias is required to strongly couple a nanomechanical resonator to a super- 
conducting qubit. However, it turns out that having such a dc bias functionality opens up 
possibilities for other heretofore difficult-to-realize quantum dynamical investigations, one 
of which we shall focus on here. 

We shall in particular investigate the quantum dynamics of the device shown in Fig. [H 
which comprises two Josephson junctions (J J) in series with a gate electrode, and where the 
source electrode to one of the JJ's contacts the center conductor of the microwave cavity, 
while the drain electrode from the other J J contacts the ground plane of the microwave cavity. 
The following section describes how the microwave cavity design allows the application of a 
dc voltage bias Vdc to the center conductor, while maintaining a very large quality factor of 
the second microwave mode to which the JJ's strongly couple ^. For not too large a Vdc 
bias, the JJ's operate in the subgap region as a "Cooper pair transistor" (CPT), where the 
dc bias generates a tunable oscillating supercurrent through the CPT via the ac Josephson 
effect. The tunneling Cooper pairs will both emit into and absorb photons from the second 
microwave mode, and it is the resulting coupled CPT-cavity mode quantum dynamics that 
will be of central interest to us. 

Related devices comprising one or more JJ's embedded in a microwave cavity date back 
to just a few years following the discovery of the ac Josephson effect where classical 
signatures of the resonant microwave modes of the tunnel junctions themselves, interacting 
with the alternating tunnel currents, were observed and discussed Q-lsl. Be ginning in the 
'90's, investigations addressed the effect of a structured electromagnetic environment with 
resonant modes on the current-voltage characteristics of dc voltage biased JJ's {q-III 



more recently, similar investigations involving double JJ devices were carried out 



. And 
12, 



131. 



However, the quality factors of the electromagnetic modes in these devices were small, typi- 



2 



cally less than 10, to be contrasted with quahty factors exceeding 10^ for the present device 
design j4| shown in Fig. [H As a consequence, emitted microwave photons will now remain 
in the cavity mode for many Cooper pair tunnel oscillation cycles before leaking out of the 
cavity mode; it does not make sense to treat the microwave cavity as an electromagnetic en- 
vironment for the CPT. Instead, the cavity and CPT should be viewed as a strongly-coupled, 
quantum coherent system. 

The CPT-cavity mode device has a mechanical analogue involving a driven coupled 
pendulum-oscillator system (Sec. IIIl B|) . The corresponding, nonlinear classical dynamical 
equations exhibit chaotic, as well as aperiodic motions depending on the initial conditions 
and the nature and strengths of the damping/noise forces. Thus, the device in principle 
allows the experimental investigation of the quantum dynamics of a system for which the 
corresponding classical dynamics is chaotic. There is a long tradition of using Josephson 
junction devices for investigating macroscopic quantum dynamics in systems with corre- 

Il5| . The Sussex group carried out some of the 



spending nonlinear classical equations jl^ . 

first, pioneering work in the '80's jiol, which was followed by the demonstration of quantum 



tunneling by the Clarke group at Berkeley [l7|, and which culminated in demonstrations 



over a decade later of superposition states by the Lukens 



18| and Mooij [l9| groups at Stony- 



brook and Delft, respectively. Subsequent, related developments have largely focussed on the 



realization of superconducting quantum bits for quantum computing applications [20|, |2l|, 
although JJ devices still occasionally are used for exploring macroscopic quantum dynamics 
and the transition to classical dynamics 22 1. 

A large body of theoretical work concerning the quantum-classical correspondence for 



driven systems has focused on the Puffing and other anharmonic oscillators 
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26 
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well as on various rigid rotor models 



40|, as 



4l|-l54| ; many insights have been gained by investigating 



dynamical properties using a quantum phase space (i.e., Wigner or Husimi function 



55 



56| ) description, and by examining the Floquet states and associated quasienergy spectra. 



However, relevant experimental results have been few 57H59|. One of the key difficulties is 
that most experimental realisations require an external ac signal to drive the system, which 
can be one of the most significant sources of noise, preventing the system from displaying 
manifest quantum dynamical behavior. In contrast, the CPT-resonator system in Fig. [1] 
generates its own ac drive, i.e., it self-oscillates. As a consequence of the ac- Josephson 
effect, only a dc voltage bias Vdc is required, and by varying Vdc, the drive frequency can 
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be tuned. Since it is considerably easier to noise filter a dc bias line than an ac-drive 
line, the device described here has considerable promise for exhibiting macroscopic quantum 
dynamical behavior. 

The outline of this chapter is as follows. In Sec. HIl we give a description of the CPT- 
cavity system. The classical system equations are derived in Sec. Illlt and solutions to 
these equations are discussed in Sec. [IVl The corresponding quantum master equation 
is derived and a quantum phase space representation of the system state is described in 
Sec. lYl Solutions to the quantum master equation are discussed in Sec. |Vl]and in Sec. IVIII 
we investigate the classical limit of the quantum master equation. We conclude in Sec. IVIIIi 

Much of the analysis will in fact deal with a simplified system comprising the driven 
'pendulum' part of the device. An analysis of the full CPT-resonator mode system dynamics, 
including results from experiment, will be published elsewhere. 

II. THE CAVITY-COOPER PAIR TRANSISTOR DEVICE 

To introduce a dc bias into a high-Q microwave cavity, we begin with a standard coplanar- 
waveguide-based resonator that is one wavelength A long at the operating frequency, illus- 
trated schematically in Fig.[T](a). As is usually the case, the ends of the cavity are terminated 
by small capacitors (on the order of a few fF) that even at a typical operating frequency of 
5 GHz have a large impedance. To a first approximation, then, we can treat these termi- 
nations as open circuits, so that the cavity voltage is a local maximum at the cavity end, 
and the cavity current a local minimum. At a distance A/4 from each end of the cavity, the 
situation is reversed: the cavity voltage is minimal and the current maximal, so that the 
A/4 points are low impedance points. 

At the A/4 low impedance points we then introduce dc bias lines consisting of sections of 
waveguide terminated with an inductance Lb. These lines are chosen to have a length A/2, 
so that the impedance they present to the main cavity line at the A/4 point is the same as 
their terminating impedance iuLf,. For even a small inductance of a few nH, this impedance 
can be substantial at the operating frequency of the cavity. 

A microwave photon approaching either dc biasing "T" junction will therefore see a short 
circuit (the low impedance of the main line) in parallel with a large impedance (the dc bias 
line) and to first order the cavity photons will be unaffected by the presence of the dc bias 
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FIG. 1: (a) Schematic diagram of a dc biased microwave cavity, showing the location of the 
inductively terminated bias lines and the sample location (black dot), (b) Illustration of a CPT 
embedded in a dc-biased cavity at the central voltage antinode. 



lines. The second (full wave) resonance of the cavity should still enjoy a very large Q of 
up to several thousand in the presence of a dc bias. By placing a CPT at the center of the 
cavity (the black dot in Fig. [Tt^a)), where there is an antinode in the cavity voltage, it should 
be possible to strongly couple the CPT to the cavity and use an applied dc bias to produce 
self oscillations of the CPT/cavity system via the ac Josephson effect. 

Electron and optical micrographs of a device based on these ideas are shown in Fig. |2l 
The cavity itself is fabricated out of a Nb film on an undoped Si substrate, as shown in 
Fig. |2t^a). Input and output lines on the left and right are coupled to the main line by small 
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FIG. 2: (a) Optical micrograph of a microwave cavity with inductively terminated dc bias lines. A 
contact pad for the CPT gate is visible at the bottom center, (b) Electron micrograph of the center 
of the cavity showing the gate lead entering the cavity at lower left, the bright Ti/Au contact pads, 
and a CPT at the right, (c) Detailed view of the CPT. The gate lead is at left the the Josephson 
junctions at the top and bottom of the central island. 



capacitors. The dc bias lines extend toward the top of the image; each is terminated by a 
small on-chip spiral inductor. Cavities based on this design have been shown to posses a 
large Q of several thousand for the full wave mode at a temperature of 4 K even when a dc 
bias voltage or current is applied to the central conductor of the cavity jj]. 

At the center of the cavity, a narrow wire to be used as a gate line for the CPT is brought 
through the ground plane of the waveguide, as shown in Fig. [2]^b). Two thin Ti/Au contact 
pads are added to the central conductor and ground plane of the cavity to the right of the 
entry point for the gate wire. These contact pads, which are driven superconducting by 
the proximity effect, allow for good metal-to-metal contact between the CPT and cavity. 
Finally, the CPT and its gate are added to the structure using standard electron beam 
lithography and shadow evaporation techniques, as in Fig. EJ^c). 
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FIG. 3: Effective circuit model of the cavity-CPT system. The cavity is modeled as a lumped 
element circuit consisting of the capacitance C and inductance L, which are in parallel with the 
CPT that is formed by two Josephson junctions in series. The circuit is controlled by the two 
voltages Vg and Vdc- Dissipation is included by incorporating an effective admittance, Y{uj), in 
parallel with C and L and an impedance Z{uj) in series with the gate voltage. 

III. CLASSICAL MODEL OF DEVICE 

A. Closed system equations 

The effective lumped element model description of the dc voltage Vdc biased microwave 
cavity-coupled Cooper pair transistor (CPT) device is illustrated in Fig. El It is supposed 
that, for the considered Vdc bias range, the CPT couples predominantly to a particular mode 
of the cavity. We neglect for the time being the cavity and CPT sources of dissipation, mod- 
eled by the parallel LC network admittance and series LC network impedance, respectively, 
focusing first on writing down the closed system equations of motion. For a typical device, 
the cavity effective capacitance C is a few pF, while the Josephson junction ( JJ) capacitance 
Cj is at least a few hundred aF, and the gate bias capacitance is about 10 aF. Furthermore, 
the effective bias line inductance Lf, is a few nH and the cavity effective inductance L is 
a few tenths of nH. Thus, the typical size hierarchies are Lf, ^ L and C ^ Cj ^ Cg. 
We shall make use of these to simplify by approximation the equations of motion. Using 
Kirchhoff 's Laws and the constitutive relations for the various lumped circuit elements, it is 
straightforward to obtain the equations of motion. In terms of the phase differences across 



7 



the two JJs, the equations are 



{C + Cj)v3i + C^2 + L'\ipi + - V^? - = 

^sm¥;i + — Lj, Vdct (1) 

and 

{C + Cj + Cg)(p2 + C(^i + L-^{}px ^^2-v\- = 

- — ^ Sin ^2 + -i^Ll^a + — (2) 

'i'O 'i'O 

where $0 = h/{2e) is the flux quantum, Ic is the JJ critical current, ip^ is an integration 
constant, and we have assumed L <^ Lf,. Transforming to 'center-of-mass' (CoM) and 
relative phase coordinates ■j± = {ipi ± '^2)/'^, Eqs. ([I]) and ([2]) become 

C7+ + i.~'(7+ - 7+) = -^sin7+cos7_ + ^L.^V^ (3) 

and 

27r/c TT 
= ^+ ~ ^'^9^9^ 

where we have assumed Cg <^ Cj C . 

Eqs. (|3]) and (jl]) follow via the Euler-Lagrange equations from the Lagrangian 

/: = 2 ^- J C7, + - ^- J C.7- + -L, Vact7. + ^C^,V.7- 

~o — ^ (7+ -7+) + — /cCOS7_cos7+. (5) 

Z V TT / TT 



The Hamiltonian is 



+Ec,(Ar-Arj^-2Ejcos7+cos7_, (6) 

where N = p^/h is minus the number of excess Cooper pairs on the island, Ng = CgVg/{2e) 
is the polarization charge induced by the applied gate voltage bias Vg in units of Cooper pair 
charge, Ecj = (2e)^/(2 ■ 2Cj) = e'^/Cj is the approximate CPT charging energy (neglecting 
Cg), i.e., the electrostatic energy cost for putting one additional Cooper pair on the CPT 
island, and Ej = 1^^0/(271) is the Josephson energy of a single J J. 

It is convenient to work instead in terms of the shifted CoM coordinate: 7+ = 7+ — 7!|. — 
cjdt, where the driving frequency is 

L eVdc .-s 



Performing this canonical transformation with the appropriate generating function, we ob- 
tain the following transformed Hamiltonian: 

H = (^y(2C)-y^+(iHy(2L)-M + i?c,(iV-iV,f 

-2Ej cos 7_ cos(7+ + Udt) , (8) 

where we have dropped the tilde on the shifted 7+ coordinate and have set 7° = 0. 

The key observation to make about Hamiltonian (j8]) is the presence of the time- dependent 
drive, which originates from the ac Josephson effect, and can be controlled via the exter- 
nally applied Vdc bias [Eq. ([7])]; the nonlinear system self-oscillates. In contrast to most 
other driven nonlinear system investigations, no externally applied ac drive is required, thus 
eliminating one of the main sources of noise that hinders the demonstration of macroscopic 
quantum dynamics. 



B. Mechanical analogue model 

In order to gain insights into the cavity mode-CPT dynamics, as well as motivate other 
parameter choices, it is useful to consider a mechanical analogue. Hamiltonian ([8]) can be 
reexpressed in the following form: 

2 

— J-W^ cos(7+ + Udt) cos 7_, (9) 

where /+ = C{<!>o/t^)^ I = h^/{2Ecj), w+ = l/^LC, and w_ = 2^EjEcj/h. From 
Eq. ([9]), we see that the cavity mode-CPT system is equivalent to a system consisting of 
two coupled rotors with moments of inertia I± and angular momentum J^. Neglecting 
the rotor coupling, the '+' rotor behaves as a torsional oscillator with frequency For 
small 7+ angular displacements and with the drive turned off (i.e., Vdc = 0), the '-' rotor 
behaves as a pendulum. For small 7_ displacements, the pendulum oscillates approximately 
harmonically with frequency a;„. With the drive turned on (i.e., Vdc 7^ 0), the pendulum 
rotor's 'gravitational acceleration' is sinusoidally modulated at frequency w^, periodically 
switching sign as a result. The gravitational acceleration is also modulated by the torsional 
oscillator coordinate. The ratio of the rotors' moments of inertia is 

1+ ^ Ri£Ecj_ ^2C_ 
9 



J, 



$0 
27r 



where Z = ix^LjC ^ 50 is the cavity impedance and Rk = h/e^ ^ 25.8 kfi is the von 
Khtzing constant. The frequency ratio for small angle, undriven oscillations is 



^ = 2 ^f^\ (11) 

For typical capacitance values, CPT charging and Josephson energies of a few Kelvins (in 
units of ks)-, and for a cavity mode frequency u;+ = 27r x 5 GHz {Tiw+jkB = 0.24 K), we see 
that the ratios in Eqs. ffTOj) and ffTTl) are large: the moment of inertia ratio is of order 10^ 
and the frequency ratio of order 10. Thus, the mechanical analogue corresponds to a fast 
pendulum with a small moment of inertia that is coupled to a slow torsional oscillator with 
a large moment of inertia. 

A measure of the zeropoint fluctuations in the pendulum angular coordinate 7_ is 



For typical CPT parameter values, we have A~p ^ 1, i.e., the zeropoint uncertainty is com- 
parable to the size of the 7_ coordinate space (= 27r radians). Thus, we don't expect the 
driven quantum pendulum dynamics to resemble much the dynamics of the driven classi- 
cal pendulum, which can be chaotic. Recovering the classical pendulum limit requires a 



60|. 



smaller charging energy than Josephson energy, for example a 'transmon'-like CPT 
The classical limit will be discussed in detail in Sec. I VIII 

How do the mechanical analogue moments of inertia compare in magnitude to those 
of actual mechanical systems? The hydrogen molecule has a rotational moment of inertia 
^ 5 X 10'^® kgm^ |61|. For ~ 5 K (= 431 /ieV = 6.9 x lO^^s j^^ ^^^^^ j_ ^ 
8 X 10~^^ kgm^. Thus, the typical CPT pendulum equivalent moment of inertia is an 
order of magnitude larger than that of the hydrogen molecule. For a transmon-like CPT, 
the moment of inertia is about two orders of magnitude larger than that of the hydrogen 
molecule. From Eq. f lTOj) . we see that the cavity mode torsional oscillator equivalent moment 
of inertia J+ is about 10^ times larger than that of the hydrogen molecule. 



C. Open system equations 



The cavity-CPT device is subject to several sources of dissipation and noise. Two signifi- 
cant electromagnetic environment sources arise from the capacitive couplings between cavity 
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and input/output microwave lines and the capacitive coupling between the CPT island and 
gate voltage bias line. Referring to Fig. [3] , we model the cavity noise/dissipation by an 
infinite parallel network of LC 'bath' oscillators, and the gate voltage noise/dissipation by 



62|-l64j . The actual dissipative mecha- 



an infinite series network of LC 'bath' oscillators 
nisms can be modeled by such infinite oscillator networks by making appropriate choices 
for the oscillator frequency distribution spectra. In the following, we will analyze the two 
noise/dissipation sources independently, beginning first with the cavity noise source. 

Extending Hamiltonian ([8]) to include the infinite parallel network of LC oscillators, we 
obtain: 

2 /.T, \ 2 



—2Ej cos 7_ cos(7+ + uat) 



2 



KI)I:|-.(|)I:^<.-..)^ 

where (pi is the phase coordinate across the network capacitance Ci. Integrating Hamilton's 
equations of motion for the network oscillator coordinate (pi, we obtain: 

=0,(O)cosa;,t + ^sina;,t + A, f d f''''''^^ ~ (14) 

where Ui = l/v^XiC^, the network oscillator "masses" are nii = C, (<l>o/(27r))' and the 
system-network oscillator couplings are Aj = 2 ($o/(2vr))^ (Lj)^^. Following the approach 



of Ref. [65|, we integrate by parts and substitute into the equations for 7+ and P4. to 
obtain the following Langevin equation: 

1 f \ ^ 2,E 

"^-^ ^ ~LC^^ " V^J ~C^ + "^'^^ ~ io "^^'^^^ ~ ^')7+(t') + fn{t), (15) 

where we have assumed that the couplings Aj are small and we have neglected frequency 
renormalization terms and where 

r(^) = (^) cosc.t = ^ E i" (16) 

is the damping kernel and 

fnit) = 7^y2^^ (m^) coswi^ + ^ sinuit] (17) 
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is the noise force. Assuming the network oscillator initial coordinates 0j(O), Pi{0) are ran- 
domly distributed according to the Maxwell-Boltzmann thermal distribution at temperature 
T, we find for the force-force correlation function: 



TT \ 1 



imUO)) = [-^J ^ksTV{t). (18) 

With ($o/7r)^C being the " mass" of the 7+ coordinate, we see that f|T8|) obeys the usual 
fiuctuation-dissipation relation. With the Markovian approximation V{t) ^ ;^'^(^)) Eq. ( ITSll 
describes a dissipative cavity mode where the admittance in Fig. |3] is simply replaced by a 
resistance R. 

Moving on now to modelling the gate voltage noise, we insert an infinite series LC network 
between the gate voltage source Vg and the gate capacitance Cg. Hamiltonian (|8]) is then 
modified approximately as follows: 



TT 



+Ecj (^N-Ng + ^Yl ~ cos 7_ cos (7+ + Udt) 



+ 

where 0j is the phase coordinate across the Ci network capacitance, and we assume Cg is 
small compared to the other capacitances. In the following, we neglect the coupling between 
the infinite series network pi and the p+ coordinates, since this results simply in adding 
to the dissipation due to the cavity mode loss considered above. Integrating Hamilton's 
equations of motion for the network oscillator coordinate Pi, we obtain: 

Pi{t) = —miUi(pi{0) sin Uit + mi(j)i{0) cos Uit 

-TTiiUiXi dt'sm[ui{t-t')]{N{t') - Ng), (20) 
Jo 

where the network oscillator frequencies and masses are the same as for the cavity mode, 
while the system- network oscillator couplings are now Aj = {2Ecj /h){Cg/Ci). Substituting 
Eq. ( 120|) into Hamilton's equation for 7_, integrating by parts and neglecting renormalization 
and shift terms, we obtain the following equation: 

7_ = ^(iV-iV,) + i Vm,A- f dt'cosuj,{t-t')N{t') 

12 



2Ecj 2E 



sin7-(^') cos(7+(t') + uj.t') + ^ Xipf\ (21) 



where 



pf\t) = —■miUJi(f)i{0) sinwjt + mi0i(O) coswjt. (22) 



Assuming the network oscillator initial coordinates 0i(O), 0i(O) are randomly distributed 
according to the Maxwell-Boltzmann thermal distribution at temperature T, we find for the 
correlation relation: 

= A:Br(^y5^^cosa;,t. (23) 
Now, we have: 

V A v^"^ _ £^ £o V- : (0) _eC^ym . 

where V^network is the fiuctuating voltage across the unloaded series network. But in the 
Markovian approximation, the voltage noise across a resistance is 

: K!e°Lk(i)K?tLk(0) ) = 2ksTRSit), (25) 

and thus 
with 

cos Uit = 2R5{t) , (27) 

where R is the effective resistance characterizing the loss associated with the gate voltage 
noise. Substituting Eq. ( 1271) into the damping term of Eq. ( 12T1) . we obtain for the 7_, 
coordinate equations in the presence of gate voltage noise and associated damping within 
the Markov approximation: 

7- = ^(iV - iV,) - ^ (^) ' i?sin7- cos(7, + u,t) + i (28) 
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and 



N = ^ sin 7_ cos (7+ + Udt). (29) 



Now that we have analyzed both the cavity noise and gate voltage noise, we finally write 
down in dimensionless form the classical Markovian Langevin equations for the cavity-CPT 
system in the presence of both noise sources. In first order form, the equations of motion 
are: 

7+ = P+ 



P+ = -7+ + / sin (7+ + UdT) cos 7_ - V+ + -^c' 



T 



7- - ^sm7_cos(7. + ...)+Ar,(.) 

IE 

N = — sin7_ cos(7+ + w^r), (30) 

where the dimensionless conversions are r = u^t and p+ = p+/(/+a;+), with w+ = 
1/VLC. The dimensionless drive force amplitude and frequency are / = ttLIc/^o = 
A{Z / Rk)Ej / {fvjjj^) and Ud = {L/ Lf,)eVdc/ {^+), respectively [with the tildes subsequently 



dropped in Eq. fl5U]) and below]. The cavity mode quality factor is Qc = Rc\/C/L in terms of 
the cavity mode resistance Rc, while Rg denotes the gate voltage resistance. The associated 
dimensionless cavity and gate bias noise "forces" satisfy the respective correlation relations 

mr)Km = 2 (^^^ LkBT,Q-'5{T) (31) 

and 

where we distinguish the cavity mode environment and gate voltage effective noise temper- 
atures, since they are not necessarily the same in experiment. 



IV. CLASSICAL DYNAMICS 



The set of Langevin equations ( l30l) provides a full description of the classical stochastic 
dynamics of the system. Numerical integration of these equations averaged over many 
different realizations of the noise allows one to obtain probability distributions for all of the 
system variables. Ultimately these distributions could then be compared with appropriately 
chosen quasiprobability distributions for the corresponding quantum degrees of freedom. 
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However, this approach is rather demanding from a computational point of view, especially 
for the quantum dynamics. We will restrict ourselves to outlining the behavior of the simpler 
system consisting of the driven Cooper-pair transistor alone. In effect this corresponds to 
the limit of small /, and Qc- 

Looking at Eq. f l30p . it is clear that for a strongly damped and weakly driven cavity, the 
variable p+ will remain small and hence to a good approximation it will be possible to drop 
the p+ dependence of the A^, 7_ equations so that the latter become entirely decoupled from 
the evolution of the cavity variables. In this limit we are left with just the pair of equations. 



2Ej 

N = —- — sin7_ cosfcjrfr). (33) 

We start by solving Eq. fl55]) in the limit where Rg = 0. In this regime the equations 
are simple classical equations of motion for and 7_. Nevertheless, they reveal a conaplex 



dynamical behavior which has already been investigated in different contexts (see e.g. |52|). 
Depending on the initial conditions and the choice of parameters, the system typically has 
a mixed phase space in which the behavior is either chaotic or quasiperiodic. The phase 
space is visualized in a stroboscopic plot in which a point is plotted after each period of the 
drive, examples of which are shown in Fig. HI In the limit Ej — )■ the system is integrable 
with natural frequencies 2Er; N, hence for very small values of Ej/Ecj the phase space is 



perturbed around resonances |53| which occur at = ±uJd/ {2Ecj) (see Fig. UK); as Ej/Ecj 
is increased the resonances get larger and a chaotic sea forms when they overlap. Islands 
of stability (where the orbits remain quasiperiodic) are found near = ±uJd/{2Ecj) even 
when Ej/Ecj > 1 (see Fig. ib). 

We can explore the sensitivity of the system to dissipation (as opposed to noise) by setting 
Tg = and changing the value of Rg. We find that even rather low levels of dissipation can 
have a significant effect on the the long time behavior. For example, for the parameters 
used in Fig. with Rg ~ 50^2, the phase space appears to contain only two attractive fixed 
points (one associated with each of the resonances). However, the chaotic sea is present as 
a transient, albeit one which can be rather long-lived: for certain initial conditions it only 
disappears after > 10'^ periods of the drive. 

Before examining the full behavior of Eq. ( 133|) with dissipation and noise, it is also worth 
considering the effect of averaging over an ensemble of initial conditions. In order to make 
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FIG. 4: Stroboscopic phase space portraits in the absence of damping. In (a) Ej = cos(1.55) w 
0.021, Ecj = sin(1.55) ?a 1.000 while in (h) Ej = cos(0.3) ?a 0.955, Eqj = sin(0.3) ^ 0.296. The 
other parameters are cod = 1, Ng = 0, Rg = 0. (For numerical calculations, energies are measured 
in units where fkvd = 1) 
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FIG. 5: Effect of a distribution of initial conditions, (a) Location in phase space of 1000 trajectories 
at T = 25r^ with initial conditions chosen from a Gaussian distribution with (N) = (7_) = 0, 
AiV = A7_ = I/V2 (b) Evolution of cos7_ for a single trajectory starting at (N) = (7-) = 0.1 
compared with an ensemble of 5000 trajectories with a Gaussian distribution of initial conditions 
centered on the same point. The parameters are Ej = cos(0.3), Ecj = sin(0.3), cod = 1, Ng = 0, 

Rg = 0. 

a comparison with the quantum dynamics we need to consider how an initial distribution of 
N, 7_ values evolves. Because of the chaotic behavior of the system the effects of considering 
a range of initial coordinates can be very dramatic even after a relatively short period of 
time. Starting from a Gaussian distribution of initial states centered on a point in the 
chaotic sea, leads to a set of trajectories that spreads out rapidly over the chaotic sea as 
can be seen in Fig. The islands within the chaotic sea stand out (the handful of points 
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FIG. 6: Stochastic probability distributions at times (a) r = Wt^i and (b) r = 25Tfi (the scale 
gives the probability of each square relative to the maximum). The parameters are Ej = cos(0.3), 
Ecj = sin(0.3), ujd = I, Ng = 0, Rg = 50017 and ksTg = 2Ej. Note we take Cg/Cj = 0.1 
throughout. 

that lie within the islands come from initial points that didn't fall within the chaotic sea). 
Clearly averaging over a range of initial conditions has a dramatic effect on the dynamics 
of the averages of the system, this is particularly clear for the quantity (cos7_) which very 
rapidly becomes a periodic oscillation with period = 2Tx/uJd as shown in Fig. [5b- 

Examples of the probability distribution for the classical, noisy, evolution of the system 
are shown in Fig. O The numerical interaction is carried out using a generalization of the 
Heun method used for deterministic differential equations 66|. In this case an average is 



carried out both over realizations of the noise and the initial conditions which are chosen 
from a Gaussian distribution with variances AA^ = A7_ = 1/v^ centered on a given 
point in phase space. When noise is added to the system the trajectories eventually diffuse 
between the chaotic sea and the quasiperiodic orbits so that the difference in the probability 
distribution over the island and chaotic sea regions gets washed out over time. In Fig. [HI 
the remnants of the island can be seen at r = lOr^, but by r = 2bTd they have disappeared 
completely. 

V. QUANTUM MODEL OF DEVICE 
A. Quantum master equation 

The Poisson bracket relations for the classical canonical coordinates are 



{7+,p+} = l; {i-,N} = h-' 
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(34) 



(where recall = p-/h). Applying the correspondence principle, the quantum commutation 
relations are 

[%,p+]=th; [^.,N]=t. (35) 

However, the phase coordinates are not periodic functions of their associated system config- 
uration spaces; the representations of the commutation relations fl35|) give the unbounded 
eigenvalue spectrum M for the corresponding phase operators. While this is not a problem 
for the 'torsional' oscillator because of the strong harmonic confining potential, which limits 
the accessible region of configuration space, the 'pendulum' typically explores the whole of 
its unit circle (S^) configuration space. A suitable pendulum configuration space function 
is e^''' with Poisson bracket relation: 

{e'^-,N} =ih-'e'^-. (36) 

The corresponding commutation relation is then 

[e^^-,iV] = -e*^-. (37) 

Eq. (1371) has infinitely many unitarily inequivalent representations jo^] that can be labelled 
by a real parameter < 5 < 1. Each representation is spanned by a number basis |A^)5, 
where 

N\N)s= iN + 5)\N)s, N = 0,±1, ±2,.... (38) 
Introduce raising and lowering operators for the torsional (CoM) coordinate: 



7+ = -l/^^^la + a^); P+ = «y ^^^^|^(a+ -a), (39) 

where recall w+ = (LC)^^/^, the 'moment of inertia' is /+ = C($o/^)^, and we have dropped 
the hats on the operators for notational convenience. The CoM phase coordinate oscillator 
zero-point uncertainty is 



where recall Z ^ 50 is the cavity impedance and Rk = /i/e^ ~ 25.8 kQ is the von Klitzing 
constant, so that = (50/25800) ~ 0.04. The Hamiltonian operator corresponding to 
dHD is 

n = hu+a+a + Ecj Yl (N + 6 - Ngf \N)s5{N\ 



N=-oo 
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-Ej {\N + I) 56{N\ + \N- I) S5{N\) COS [^%{a + a-^)+ujdt\. 

(41) 



A''=-oo 



The parameter 5 appearing in the Hamiltonian operator is a purely quantum signature of the 
nontrivial topology of the corresponding classical pendulum's configuration space 5"^. An 



interesting question concerns the particular value for 6 that Nature chooses and why [67h69|. 
However, it is likely not possible to measure 5 in experiment, since from (HTl) it is clear that 
the effect of a nonzero 5 value is indistinguishable from that due to the presence of an excess 
charge on the CPT island. From now on, we shall set 5 = 0. 

We now derive the open system quantum master equations within the self-consistent Born 



approximation (SCBA) following the approach reviewed in Ref. [70|. In the following, we 
analyze the two noise/dissipation sources independently, beginning first with the cavity mode 
environment. We write the Hamiltonian ( fT3|) as "H = l-Ls + 'He + y where the Hamiltonian 
Tis describes the cavity mode-CPT system [Eq. fj4T]) ]. the Hamiltonian T-Le describes the 
infinite parallel LC network environment, and the interaction part is 

^ ^ i * i 

where Aj = 2($o/(27r))^(Li)^"'^. Defining B = we obtain for the master equation 

within the SCBA: 

pit) = -{[V-s.p{t)] 

1 n 



+ \{[B{t\B{t')]) [7+, {7+(t' - t),p(t)}]} , (43) 

where the operators B{t) and 7+(t' — t) are in the interaction picture and the expectation 
values are performed assuming the environment (infinite parallel LC network) is in a thermal 
state. We have 

-({fi(t),5(0)}) = - / duj. J {uj) cos ujt (I + 2N{u)) 
2 TT Jo 

1 i 

-([5(t),5(0)]) = / dojJ{uj)smojt, (44) 

2 vr Jo 

where N[uj) = [exp{hw/ ksTc) — is the thermal occupation number of the environment 
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at frequency u and where the spectral function is 



with rrii = Cj(<l>o/(27r))2 and cOi = 1/a/L^. 

In Sec. IVIIt we compare the quantum versus classical dynamics and establish conditions 
under which the former is well approximated by the latter-the so-called classical limit. A 
necessary condition to be in the classical limit is that the environment temperature must be 
sufficiently large such that we can make the approximation A^(ci;) ~ ksTi./ {fvjj) 3> 1. This 
requires ksTc ^ /la;+. The environment correlation function then becomes 

h{B{t),Bm)=kBT,y^^co^uj,t=(^\ CkBT,T{t), (46) 

where r{t) is the classical damping kernel f|T6|) . If, furthermore, the spectral function upper 
cut-off satisfies A ^ ksT^, then we can make the Markovian approximation T(t) ^ '^^i'^) 
and 

Substituting expressions (H6ll and ( H71) into Eq. (l43l) . integrating by parts and using p+ = 
($o/7r)^ (^7+, we obtain 

m = [Hs,p{t)] - [7+, K,p(t)}] - y^I+ksT, [7+ [7+,p(t)]] , (48) 

where F = 1/{RC) is the damping rate. Eq. ( HHll is just the standard Born-Markov master 
equation for a quantum Brownian particle in the high temperature limit [7o|, where the 
second term on the right hand side describes damping and the third term on the right hand 
side describes diffusion. 

While Eq. f HHj) is appropriate for investigating the classical limit, under the cryogenic 
conditions of an actual experiment and for say an a;+ > 27r x 5 GHz cavity mode, we expect 
that ksTc ^ fvjJ^, so that a low temperature limit is more appropriate. Using Eq. ( I39l) to 
express the master equation in terms of raising and lowering operators, making the rotating 
wave approximation and the replacement /cbTc — > /ia;+/2, we obtain the following 'low 
temperature' master equation: 

i 1 

p = —t[^s,p] — (a"'"ap + pa'^a — 2apa'^) . (49) 
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The non-Hermitian part of the master equation (H9!) is of the Lindblad form, ensuring that 
the solution to (1491) for the density matrix p{t) is consistent. However, in the above outhned 
derivation of this master equation, there is no apparent justification for the replacement 
fc^Tc — 7- since is not in general the characteristic frequency for dynamics de- 

scribed by the nonlinear, time-dependent Hamiltonian Tis- Only for a harmonic oscillator 
Hamiltonian with classical frequency is this replacement justified, provided the damping 
is sufficiently weak, i.e., KT <^ ksTc ^ hwj^. Nevertheless, for better or worse, we shall 
follow common practise and assume that master equation (H9!) provides an adequate model 
for weakly damped cavity mode-CPT quantum dynamics at low temperatures. 

We move on now to derive the master equation modelling gate voltage noise, starting from 
the Hamiltonian flTI?]) . For the master equation within the SCBA, we obtain [c.f., Eq. fH5]) ] 

m = ~[ns,pit)] 

+ \mt).B{t')]) [iV,{iV(t' -t),p(t)}]} , (50) 

where now B = Yli^iVii with Aj = (2Ecj/h){Cg/Ci). The correlation functions still take 
the same form as Eq. (144|) . but where now the spectral function is 

= (51) 

i 

In the high temperature limit, appropriate for comparing the quantum versus classical dy- 
namics, the environment correlation functions become 

2 T / ^ \ 2 



l{{B{t),Bm) = fcBT,(^y5^^cosa;,t = 2A:^T(^yi?,5(t) 



i({i?(t),i?(0)}) 



(52) 



where the last equality results from making the Markovian approximation. Substituting the 
correlation relations (!52|) into the master equation ( |50|) . integrating the damping term by 
parts and using also Heisenberg's equation to solve for A^, we obtain: 

m = ~ [Us, p{t)] - ^ ' ksT.R, [iV, [iV, p{t)]] 
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Eji?g[Ar,{sin7_cos(7+ + a;dt),p(t)}], (53) 

where the second term on the right hand side describes diffusion and the third term describes 
damping. Note the atypical, exphcit time-dependence in the damping term. 

In contrast with the more famihar quantum Brownian master equation ( l48l) for the cavity 
mode environment, there is no corresponding simple prescription for recovering from Eq. ( 153|1 
a consistent, low temperature master equation that is appropriate for lower gate voltage noise 
levels expected in an actual cryogenic experiment, where the bias lines are filtered. Given 
the difficulties in finding such a low temperature master equation, we shall instead in Sec. I VII 
'take the path of least resistance' and simply compare the quantum dynamics in both the 
presence and absence of the environment non-Hermitian terms in ( l53l) . so as to gain some 
understanding of the system quantum dynamics in the presence of gate voltage noise. 




B. Wigner functions and coherent states 

In order to make a comparison between the quantum and classical dynamics of our sys- 
tems we need to identify suitable tools with which to describe the quantum dynamics. For 
continuous systems coherent states prove useful as initial states because they are localized in 
phase space and have minimum uncertainty making them the quantum states most closely 
connected to a classical phase space point. Furthermore, the Wigner function quasiprobabil- 
ity distribution provides an effective way of visualizing the quantum evolution of a continuous 
variable system in phase space and signals the presence of quantum interference effects by 
turning negative. Whilst both conventional coherent states and the Wigner function can 
be applied directly to the cavity degrees of freedom, the Cooper-pair transistor is different 
since the Cooper pair number is discrete. Nevertheless one can define appropriate versions 
of both coherent states and Wigner functions for the N, 7_ degrees of freedom, but with 
some important differences compared to the usual continuous variable case. 

Given the basic algebra of our number and phase operators, [7_, A^] = i and [A^, e'^-] = 



e*'''", we can use a definition of coherent states first deve^ 
rotation angle variables. We define the coherent states 68 



oped for angular momentum and 



69j as eigenstates of the operator 
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X = e*^" ^ and hence they take the form 



^ +00 
1^) = n^Yl 



j=-oo 



where the normahzation factor is 



+00 

j=-oo 



The complex parameter x can be written in terms of a charge N and phase 7_, x = 
which are closely related to the expectation values of the corresponding operators 



68| 



{X\N\X) ^ N 



and 



(X|e^^-|x) 



(54) 



(55) 



Although the properties of these states differ somewhat from those of the harmonic os- 
cillator coherent states they nevertheless take a fairly simple and intuitive form in phase 
space. Again, following the approach developed for angular momentum and rotation angle 



variables, we can use the form of the Wigner function developed for this case 



71 



73|. For a 



system with density operator p the quasiprobability distribution is given by [73 1 



IT 



(7- + 7-|p|7- -7-) 



-7r/2 



sin[(iV- A^'-/i/2)7r] 



^=0,1 Ar'=-oo ^ ^' ' 



WN'+^,/2h-,t), 



(56) 



where 



WN+fj./2 



-, +00 

(7_,t) = - y e-''^'''-^''~/'^''~{N-N'\p\N + N' + ij,) 

TT f ^ 



N'=-oo 



^^^e-2^(7v+M-/2)7L (7_ - 7^|p|7_ + 7^) 



(57) 



and |7_ + 7^) is one of the phase states. Note that the w functions have definite parity, 
i.e., WN+f_i/2{l- + 7r,t) = (— l)^WAr+^/2(7-, t). Thus the Wigner function takes the form of a 
discrete series of strips labeled by the Cooper-pair number and which are continuous in the 
phase. 
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FIG. 7: Wigner functions for a charge state = 3) (a) and a coherent states \x = e*'^/^^^) (b). 
Note that m (b) there are regions around 7_ = —tt/2 and = 1,3 where the Wigner function 
becomes very shghtly negative. 

Examples of Wigner functions for a charge state = 3) and the coherent state \x = 
gi7r/2-2^ are shown in Fig. [71 The Wigner function for a charge state, = M), takes a 
very simple form: W^Ar(7_) = 5Af,j\//(27r), i.e., it is just a uniform strip for = M and zero 
elsewhere. Interestingly the charge states are the only examples of pure states which give 
rise to Wigner functions which are positive everywhere {tJ] . For number-phase variables the 
Wigner functions of coherent states do have negative regions, in contrast with the situation 
for continuous variable systems However, if we wish to consider an initial quantum 

state which is analogous to a classical point in phase space the coherent states are still a 
very good choice as the amount of negativity which actually occurs is in fact very small in 
practice [See Fig. [7b] and our only choice if we wished to eliminate the negativity entirely 
while still using a pure state is to work with charge states which are completely spread out 
in phase. As we can see from Fig. [7)d the Wigner function for \x = e*'^/^^^) is strongly 
peaked around = 2,7_ = 7i/2 and, apart from very small regions of negativity [3], is 
very reminiscent of the corresponding continuous variable case. 



VI. QUANTUM DYNAMICS 

In exploring the quantum dynamics we again choose to focus on just the behavior of 
the Cooper-pair transistor, as we did with the classical dynamics. This in effect means 
that we take the limit — )■ in the Hamiltonian [Eq. ( 14T]) ]. The presence in the model 
of a non-zero gate resistance and large effective gate voltage noise temperature ensures 
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that the long-time behavior of the CPT system is classical for sufficiently large Ej, in the 
sense that the Wigner function will be everywhere positive and also a smoothly varying 
function of (see Sec. IVIip . However, over short times a very different picture emerges: 
even very 'classical' choices of the initial state can evolve naturally into states with strongly 
non-classical features. 

We start by considering how the system evolves starting from an initial state which is 
'classical' in the sense that its Wigner function is relatively smooth as a function of A^. The 
initial state we choose to use (illustrated in Fig. [8^) is the steady-state of the system when 
the dc voltage is switched off (i.e. the Ud = limit), and hence should be easy to prepare 
in practice. Figure [H] shows a series of snapshots of the Wigner function at different times 
after we set = 1 at time t = 0. The Wigner function becomes stretched, reaches the edge 
of the phase space, and then starts to warp around on itself. This wrapping around leads 
to interference and the formation of regions where the Wigner function is negative. Very 
similar results are obtained when the system is initially in a coherent state instead. 

Over longer time-scales and for sufficiently large Ej / Ecj the quantum dynamics matches 
up fairly well with the stochastic dynamics (described in Sec. IIVI) . A good overview of the 
dynamics can be obtained by looking at some of the moments as a function of time, as shown 
in Fig. [9l The behavior of (A^^) and (cos7_) in both the quantum and stochastic dynamics 
is eventually periodic, a consequence of the underlying periodicity of the drive. The average 
charge, (A^) decays rapidly to zero in both cases. The classical stochastic dynamics comes 
close to reproducing the behavior seen in the quantum dynamics of both (A^^) and (cos7_) 
in the long time limit, although the amplitude of the stochastic oscillations is slightly smaller 
than the quantum ones. Going beyond the moments of the system, we can compare the full 
stochastic probability distribution with the Wigner function of the corresponding quantum 
evolution as shown in Fig. [TOl 

Different choices of initial conditions corresponding to different regions of the classi- 
cal phase space lead to rather different quantum evolutions, although again this is only 



a transient effect in the presence of dissipation. Analyses of similar systems 5ll, l52 | 



(wit 
ing 



lout dissipation) 



27 



48 



51 



52 



57 



lave shown that they can display chaos-assisted quantum tunnel- 
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77 



79| . In the classical dynamics any point initially chosen to be 



within one of the stable islands is confined to either the A^ > or A^ < part of the phase 
space (depending on the initial conditions) and the value of A^ oscillates quasi-periodically 
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FIG. 8: Spreading and interference during the early stages of quantum evolution. Wigner functions 
for a system where ujd is switched from to 1 at t = 0: (a) r/rd = 0, (b) r/r^ = 0.477, (c) 
r/Td = 0.716, (d) r/Td = 0.955. The gate resistance is 50f], with other parameters Ej = cos(0.2), 
Ecj = sin(0.2), Ng = 0, and k^Tg = 2Ej. The initial state is the steady-state of the system with 
ujd = 0. 




FIG. 9: Comparison of the dynamics of the stochastic and quantum evolution (a) (A^^) (b) 
(cos(7_)). The parameters chosen are: Ej = 4, Ecj = 1, Ng = 0, ksTg = 2Ej and Rg = 5000. 
The initial state for the quantum (stochastic) dynamics was a coherent state (Gaussian distribu- 
tion) centered on the point = 7_ = 1. 
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FIG. 10: (a) Stochastic probability distribution and (b) Wigner distribution of the corresponding 
quantum dynamics at r = 75rrf. The parameters are the same as in Fig. [9l 

on a timescale In the quantum dynamics, the system can tunnel between these two 

regions of phase space. We looked at the evolution of an initial coherent state centered on 
points which are either within a stable island or the chaotic sea in the classical phase space 
(with parameters which match those of Fig. |Dd). The dynamics of the average charge for 
these two cases is compared in Fig. [11] An initial state centred on a stable island (Fig. 
[TTk) shows fast oscillations in (A^), which have the same time-scale and (initially) a similar 
amplitude to those seen in the classical dynamics. However, there is also an underlying 
much slower oscillation which takes the system to the opposite side of the phase space. In 
contrast, for an initial coherent state centered in the chaotic sea (Fig. ITTb) the behavior is 
much less regular. The differences in behavior wash out very rapidly when dissipation is 
included, an important sign that the slow oscillations involve coherent superpositions. The 
slow oscillations in Fig. [TTb take the system from a state localized around one of the regions 
of phase space corresponding to a classical island of stability, to one which is largely localized 
on the other stable island. The system tunnels via a state which shows strong interference 
effects, as can be seen from the Wigner functions calculated after a quarter and after a half 
of the slow oscillation period shown in Fig. [T21 

VII. CLASSICAL LIMIT 

A Wigner function's domain is given by the eigenvalue spectra of the canonically conjugate 
operators from which the function is constructed. For the number/phase Wigner function 
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FIG. 11: Comparison of the evolution of (N) with (Rg = 5017) and without damping for an initial 
coherent state centered on a point (a) inside an island of the classical phase space {N, 7_) = (2.27, 0) 
(b) inside the chaotic sea (N,^^) = (2.27, 1.25). The parameters used match those in Fig. (Da. 




FIG. 12: Wigner functions at (a) r = 50t^ and (b) r = lOOr^^ for the parameters corresponding to 
Fig. dBi with Rg = 0. 

(see Sec. IV Bl) . the domain consists of parallel lines [7- = (— vr, tt) and = 0, ±1, ±2, . . .]. 
Furthermore, since these eigenvalues are the possible outcomes of measurements of the as- 
sociated observables, it is natural to use the Wigner function to define a quantum phase 
space dynamics from which the classical phase space limiting dynamics can be recovered 
by approximation. In the following, we again restrict ourselves to the driven pendulum 
subsystem dynamics only, resulting from formally setting = in the Hamiltonian (HI]). 

In contrast to a nonlinear system with the configuration space topology M, such as the 
commonly investigated one-dimensional Duffing oscillator, it does not appear possible to 
derive a closed form equation for the Wigner function W^i'^-yt) starting from the master 
equation (|53|) . The problem lies in the sin7_ potential and damping terms, which give rise 
formally to Wigner functions shifted by 1/2 in their A^- arguments. Because of the non-trivial 
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configuration space topology, the coordinate domain of tlie Wigner function consists 
of the discrete integers, not half-integers. On the other hand, it is possible to write down 
closed form equations for the definite parity functions WAr+^/2(7-, t). From ( 13^ and 
we obtain after some algebra: 



dt h 97_ 

sin- 



2,E 

■-^ sin7_ cos{uJdt) {wN-l/2+^,/2 - 



^9 



1 feCnV „ „ d 



^2 ^J^9q;:^ [sin7_ cos(a;rft) {wN-i/2+ti/2 + W7V+1/2+/./2)] ■ 

(58) 

Using the 7_ integral form of the Wigner function in ( l56l) to analytically continue the 
function to arbitrary real N, we obtain from Eqs. fl58l) and fl56|l the following equation for 
the Wigner function: 



2Ej 

h 



sin7_ cos(wdi) - W^iv+1/2) 



^9 



1 /ea\^ „ „ 9 



f^2 \Q I ^^^3^^ [sin7_ cos(a;dt) (M/Ar-i/2 + WAr+1/2)] • 

(59) 

Under conditions where W varies slowly with A^, we can Taylor expand the fractionally- 
shifted W functions to first order as a good approximation and obtain the following classical 
master equation for the probability distribution P(7_, A^, t): 



dP 2Ecr,,, , dP 2Ej . , ,dP 



2 



2 / eC \ d 
-^[-cfj EjR,^[sm^.cosiu,t)P]. (60) 

This master equation is equivalent to the classical pendulum Langevin equation fl30l) . with 
(7+,p+) set to zero. 
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In contrast to the usual situation for a nonlinear system with trivial configuration space 
topology M, the classical limit f l60|) of the quantum master equations fl58|) and f l59|) was not 
obtained by identifying and then discarding a higher derivative quantum term involving 
the anharmonic system potential [24]. Rather, the difference between the quantum and 
classical equations is more subtle and linked to the discreteness of the number (equivalently 
angular momentum) operator, which in turn is a consequence of the non-trivial configuration 
space topology. All that is required to recover the classical master equation to a good 
approximation is that the Wigner function varies by only a small amount as its argument 
N increases or decreases by one. 

Clearly, the pendulum state must have non-negligible overlap with a large number of 
angular momentum eigenstates |A^) if the Wigner function is to depend smoothly on A^; 
recovering the pendulum classical limit necessarily requires Ecj ^ Ej, e.g., a 'transmon'- 



like CPT 



60|. However, the latter inequality is not a sufficient condition: as shown in 



Fig. IH], an initially 'classical', i.e., smooth, positive practically everywhere Wigner function 
can evolve through stretching and shrinking into a Wigner function that is non-smooth in N, 
so that the classical master equation approximation ( 160|) breaks down. Furthermore, because 
of the periodicity of the configuration space, the stretching pendulum wavefunction can 
eventually interfere with itself, resulting in an oscillatory Wigner function with significant 
negative regions. 

Even though the initial state-dependent 'transient' dynamics will display quantum fea- 
tures, the dynamics will eventually settle into a steady state where the Wigner function is 
practically positive everywhere and well-approximated by the classical master equation (I6OII . 
provided ksTg >^ Ej ^ ^Cj- The smaller is the gate voltage resistance Rg (or gate capaci- 
tance Cg), the longer is the duration of the transient quantum interval. 

Interestingly, when the charging energy is not small, i.e., Ec, > -Ej, the steady state 
Wigner distribution will still be practically positive everywhere and hence interpretable as 
a probability density, provided ksTg > Ecj- However, because the Wigner function is 
non-smooth and with non-negligible support over only a small range in A^, the classical 
pendulum master equation no longer accurately describes the Wigner function dynamics. 
The question then arises as to whether there is an approximate classical description that 
is distinct from the classical pendulum equation. Such an equation must necessarily treat 
A^ as a discrete coordinate and so is more appropriately interpreted in terms of the charge 
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dynamics of the CPT. For Ej <^ Ec, and sufficiently large gate voltage resistance and 
effective noise temperature, P{E) theory may provide an adequate classical description in 



terms of incoherently tunneling Cooper pairs 



10 



80 



81| . Otherwise, Cooper pair tunneling 



across the Josephson junction is an inherently quantum coherent process: even though 
the steady state Wigner function evolution is practically positive everywhere, the model 
dynamics must be interpreted as quantum in nature when the charging energy is large. 



VIII. CONCLUSION 



In this chapter, we have investigated a strongly non-linear superconducting device con- 
sisting of a Cooper pair transistor (CPT) that is coupled to a dc voltage biased microwave 
cavity and driven by the dc bias via the ac Josephson effect. Our main focus has been on 
comparing the quantum and classical dynamics of the system - the " quantum-classical cor- 
respondence" - in particular establishing the circumstances under which the corresponding 
dynamics are similar. 

We derived the corresponding classical Langevin and quantum master equations of mo- 
tion, which describe the dynamics of the CPT-cavity system in the presence of an envi- 
ronment consisting of dissipative circuit elements. Although we did not investigate the 
dynamics of the full system, we did explore the dynamics of the driven CPT (which is the 
non-linear element of the device) in the limit where the cavity does not affect its behavior. 
The quantum-classical correspondence was elucidated by considering the Wigner function 
representation of the evolving CPT quantum state. Numerical simulation of the quantum 
and classical dynamics shows that, apart from initial state-dependent transients, the evolu- 
tion becomes very similar in the limit of large Josephson energies, for which the discreteness 
of the Cooper pair number on the CPT island is unimportant. Interestingly, though, the 
transient behavior of the system can lead to highly non-classical states even when the initial 
state is apparently very classical. 

Future work will need to explore how the full system (CPT and cavity degrees of freedom) 
behaves, as well as connect the dynamics of the system to quantities that are measured in 
experiment. On a more technical level, further analysis of the interaction between the CPT 
degrees of freedom and the gate impedance is needed in order to derive the correct description 
of the dissipative dynamics in the low temperature limit. 
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